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1.  Introduction 


Water  vapor  in  the  atmosphere  may  be  considered  as  originating  from  two  sources:  evaporation 
from  the  water  surface  associated  with  the  volumetric  water  content  of  the  soil  and  the  molecular 
diffusion  of  vapor  from  the  remainder  of  the  soil  pores.  The  water  vapor  content  of  the 
atmosphere  must  be  considered  when  dealing  with  atmospheric  processes  involving  the  turbulent 
kinetic  energy  budget  or  the  surface  energy  balance.  Specifically,  when  dealing  with  buoyant 
heat  fluxes  in  the  surface  boundary  layer,  the  portion  of  the  energy  contributed  by  the  latent  heat 
flux  must  be  accurately  determined;  otherwise,  uncertainties  will  occur  in  subsequent 
calculations. 

Evaporation  from  a  bare  surface  or  evapotranspiration  from  a  vegetated  surface  is  a  function  of 
soil  types,  soil  heat  fluxes,  net  radiation,  volumetric  water  content  of  the  soil,  and  the  time  since 
natural  or  artificial  irrigation  ended.  Natural  irrigation  may  be  defined  as  precipitation. 
Parameterization  of  the  static  or  dynamic  structure  of  the  surface  boundary  layer  of  the 
atmosphere  can  be  accomplished  by  using  semiempirical  or  engineering  estimate  approaches. 
These  techniques  can  be  based  upon  much  simplified  and  modified  forms  of  the  Penman  (1948) 
potential  evaporation  equation,  which  can  be  combined  with  the  dynamic  similarity  of  flows 
theory  of  Obukhov  (1946). 

2.  Background 

Semiempirical  parameterization  and  amplification  of  any  scheme  leading  to  practical  engineering 
estimates  of  atmospheric  surface  boundary  layer  processes  require  some  clarification  and  the 
tying  together  of  some  seemingly  unrelated  principles.  A  primary  assumption,  based  upon  the 
work  of  Brunner  (1977),  is  that  the  relative  humidity  is  constant  with  height  in  the  first  100  or 
so  meters  of  the  atmosphere.  This  premise  allows  specific  humidity  and  vapor  pressure  profiles 
to  be  established  from  wet-  and  dry-bulb  temperature  measurements  at  one  level  above  the 
surface. 

The  saturation  vapor  pressure  e,  can  be  calculated  for  any  temperature  T  at  any  height  z  using 
Teten’s  (1930)  formula  in  the  modified  form 


eg  =  6.11  exp 


17.4(7*  -  273.16) 
T  -  34.16 


(1) 


where  T  is  in  degrees  kelvin. 
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If  the  relative  humidity  is  known,  then  the  ambient  vapor  pressure  is  found  from 


e-e(  (RH) 


(2) 


where  RH  is  the  relative  humidity.  Specific  humidities  are  determined  from 


4 


0.621 


a 

p  -  0.379a 


(3) 


where  p  is  atmospheric  pressure  in  millibars.  The  density  of  moist  air  may  be  found  from 


Pw 


(4) 


where  R  is  the  gas  constant  and  0„  the  virtual  potential  temperature  given  by 


0. 


0 

1  -  0.379e/P* 


(5) 


Actual  evaporation  rates  influence  not  only  the  latent  heat  flux  into  the  atmosphere,  but  also 
buoyancy,  dynamic  stability  and,  of  course,  the  specific  humidity  profile,  which  is  typically 
written  as 


Si 

* 


4  -4, 


In  ±  (|) 


(6) 
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when  q*  is  a  scaling  humidity,  k  is  Karman’s  constant,  q  specific  humidity,  q„  the  specific 
humidity  at  z,  the  surface  roughness  length,  z  height  and  1fB  (zJL)  the  diabatic  influence  function 
for  heat.  In  the  Obukhov  (1946)  similarity  theory,  heat  and  mass  transfer  are  taken  to  be 
identical,  hence  the  use  of  ifg  (zJL)  for  moisture  profiles.  According  to  Myrup’s  (1969) 
semiempirical  results, 


RH 

1000 


3.74  +  2.67 


where  T0  is  the  temperature  at  z,  in  degrees  Celsius. 

The  surface  energy  balance  is  usually  written  in  the  form 


rn  =  h  +  le  +  g 


(7) 


(8) 


where  R*  is  the  net  radiation,  H  the  sensible  heat  flux,  LE  the  latent  heat  flux  and  G  the  soil 
heat  flux.  The  sensible  heat  flux  with  respect  to  the  eddy  diffiisivity  is  given  by 

H  *  -cp  Kg  (9) 

dz 

where  Cp  is  the  specific  heat  of  air  at  constant  pressure  (c,,  =  1014  J  k,1  “K1,  p  density,  KH 
the  eddy  diffusivity  and  0  potential  temperature.  Equation  (9)  can  also  be  written  in  terms  of 
U.  and  T*  as 


H  =  ~cp  put  r 


GO) 


or 


H  = 


-Cppw  u,  0* 


(ID 
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where  u.  is  the  friction  velocity,  T*  a  scaling  temperature  and  0*  a  virtual  scaling  temperature. 
Similarly,  the  latent  heat  flux  (LE)  is  given  by 

L  E  «  -  p  L  n.  q0  (12) 


where  L  is  the  latent  heat  of  vaporization  and  E  the  evaporation  rate.  The  buoyant  heat  flux  is 
now  written  as 


H'  •(H  +  0.07  L  E) 


(13) 


3.  Evaporation 

The  release  of  water  vapor  from  the  soil  to  the  atmosphere  can  be  thought  of  as  bare  surface 
evaporation  and  evapotranspiration  from  vegetated  surface  and  is  a  most  complex  boundary  layer 
phenomenon.  According  to  Penman  (1948)  the  potential  latent  heat  flux  can  be  shown  to  be 


LEf- 


s  +  a 


s  +  a 


(14) 


where  S  is  the  slope  of  the  saturation  specific  humidity  curve  versus  temperature,  o  is  the 
psychrometric  constant,  e,  is  the  saturation  vapor  pressure  for  ambient  air,  ed  is  the  saturation 
vapor  pressure  at  the  dew  point,  and  /  (v^  is  a  function  of  the  horizontal  mean  wind  speed  v , 
given  by 


/®  =  °.35(1+1L) 


(15) 


over  a  grassy  surface. 
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To  establish  the  actual  evaporation  a  moisture  availability  coefficient  M  such  as  suggested  by 
Nappo  (1975)  may  be  written  as  the  ratio  of  actual  evaporation  to  potential  evaporation  given 
by 


M  -EfE,. 


(16) 


Jackson,  Idso,  and  Reginato  (1976)  were  able  to  express  M  in  terms  of  the  increasing  albedo  of 
drying  soils,  as 


M  »  («rf  -  a)l(a4  -  aj 


(17) 


where  ad  is  the  dry  soil  albedo,  aw  is  wet  soil  albedo,  and  a  is  the  albedo  for  a  particular 
time.  Albedo  measurements  were  also  found  to  be  an  ideal  mechanism  of  integrating  various 
drying  ratios  and  partitioning  the  portions  contributing  to  the  energy  limiting,  or  potential 
evaporation,  and  the  soil  limiting  phases  of  the  process.  The  coefficient  M  will  vary  from  unity 
for  a  soil  at  field  capacity  to  zero  for  a  dry  surface.  The  potential  evaporation  rate  can  be 
combined  with  the  soil  limiting  rate  E,  to  obtain  the  actual  evaporation  ratio. 

E  =  M  Ep  +  ( 1  -  At)  Et  (1*) 


Ritchie  (1972)  found  the  soil  limiting  rate  to  be 


es  -  c  rm 


(19) 


where  t  is  the  time  in  days  and  c  varies  with  soil  type  and  season  of  the  year.  Seasonal 
dependencies  are  a  function  of  temperature  and  will  vary  by  about  a  factor  of  2  from  winter  to 
summer  as  shown  in  figure  1. 

The  complexity  of  the  evaporation  process  requires  that  the  energy  limiting  and  soil  limiting 
phases  be  partitioned  according  to  the  starting  times  of  the  soil  limiting  fraction.  Thus,  equation 
(18)  must  be  rewritten  as 


E  =  MEp  +  C  Y,  (Aftl  -  Afj)  (»  -  i  +  1  rw  (20) 
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where  n  is  the  number  of  days  after  evaporation  began.  Multiplying  both  sides  of  equation  (20) 
by  L  yields  the  latent  heat  flux. 


The  moisture  availability  coefficient  M  has  been  evaluated  using  evaporation  data  extracted  from 
Dugdale  (1989),  Ritchie  (1972),  and  Jackson  (1973).  Potential  evaporation  Ep  was  taken  to  be 
the  evaporation  rate  on  day  1  after  irrigation  or  rainfall  ended  and  used  as  the  numerator  in 
equation  (16)  for  normalizing  the  data.  The  three  data  samples  were  from  widely  different 
climatic  regions  (Southern  Arizona,  Central  Texas,  and  the  African  Sahel)  and  it  is  noteworthy 
that  when  normalized  a  single  curve  fits  all  the  data  quite  well,  as  shown  in  figure  2. 

Dugdale’s  data  were  for  three  separate  rainfall  amounts,  3,  7.5,  and  15  mm;  while  the  Jackson 
et  al.  (1976)  samples  were  for  all  four  seasons  of  the  year.  Ritchie’s  data  were  for  four  soil 
types  (Adalanto  Clay  Loam,  Yolo  Loam,  Houston  Black  Clay  and  Plainfield  sand).  These 
results  although  not  definitive  suggest  that,  discounting  extremes  such  as  Jackson’s  summertime 
evaporation  rates  and  the  one  very  light  rainfall  (3  mm)  reported  by  Dugdale,  the  coefficient  M 
may  well  be  independent  of  season,  rainfall  amounts,  and  soil  type. 

4.  Discussion 

The  importance  of  considering  evaporative  processes  and  the  evaluation  of  specific  humidity 
profiles  in  the  surface  boundary  layer  significantly  concerns  water  vapor  effects  upon  the 
dynamic  stability  of  the  atmosphere  and  the  surface  energy  balance.  Atmospheric  stability  in 
the  surface  boundary  layer  is  generally  designated  by  the  Obukhov  (1946)  length,  and  according 
to  Busch  (1973)  is  given  by 


L  = 

kge* 


where  g  is  the  gravitational  acceleration. 

Richardson  (1920)  and  Obukhov  initially  ignored  the  contribution  of  latent  heat  to  buoyancy,  and 
consequently,  the  accepted  similarity  definitions  of 


z  m  _  kgzH 


L 


(22) 


where  4>„  and  4>m  are  dimensionless  lapse  rates  and  wind  shears  respectively  and  is  the  eddy 
viscosity  (only  valid  in  totally  dry  air). 
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Figure  3  shows  a  comparison  of  equations  (21)  and  (22)  in  terms  of  the  scaling  ratio  z/L.  The 
data  utilized  to  prepare  figure  3  was  extracted  from  a  compilation  by  Barad  (1958)  of  wind, 
temperature,  and  vapor  pressure  profiles  observed  during  Project  Prairie  Grass.  These  data  woe 
carefully  screened  and  12  reasonably  stationary  and  homogeneous  profiles  were  selected  for 
analysis.  Note  in  figure  3  that  the  exclusion  of  water  vapor  in  equation  (22)  indicates  an 
atmosphere  less  stable  than  it  reaily  is. 

The  surface  energy  balance,  as  given  by  equation  (8),  is  also  dramatically  modified  by  an 
increase  in  atmospheric  water  vapor.  Two  simulated  cases  are  given  in  table  1  for  relative 
humidities  of  25  and  75  percent.  Initial  conditions  are  presumed  to  be  identical  with  the 
exception  of  the  relative  humidity.  Tripling  the  relative  humidity  drastically  alters  the  surface 
energy  balance  as  a  function  of  incoming  solar  radiation  reaching  the  earth’s  surface.  Insolation 
is  reduced  by  33  W  nr2  or  about  4  percent,  while  the  net  radiation  increases  by  approximately 
4  percent.  The  most  significant  modification  to  the  surface  boundary  layer  occurs  in  the  stability 
parameters  L  and  z/L,  both  of  which  decrease  by  37  percent.  This  decrease  is  attributed  to  the 
approximately  17  percent  increase  in  the  latent  heat  flux  and  an  extremely  large,  43  percent, 
decrease  in  buoyant  heat  flux.  The  results  show  that  atmospheric  water  vapor,  or 
evapotranspiration,  cannot  be  ignored  in  the  surface  boundary  layer. 

The  increase  of  water  vapor  in  the  atmosphere  resulting  from  the  increase  in  the  relative 
humidity  has  the  effect  of  reducing  the  transparency  to  radiation,  especially  the  visible  radiation. 
The  presence  of  water  vapor  as  a  suspensoid  will  account  for  the  indicated  changes  in  the  surface 
energy  balance. 
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Table  1.  The  Surface  Energy  Balance  and  Atmospheric  Stability  as  a 
Function  of  Relative  Humidity 


Parameter 

RH  =  25% 

RH  =  75% 

Height 

2  m 

2  m 

Roughness  length 

0.0065  m 

0.0065  m 

Temperature 

27  °C 

27  °C 

Wind  speed 

4  m  s'1 

4  m  s'1 

Insolation 

875  W  m*2 

842  W  m-2 

I  Net  radiation 

675  W  nr2 

740  W  m*2 

|  Soil  heat  flux 

66Wm2 

141  W  m'2 

Sensible  heat  flux 

249  W  nr2 

127  W  m*2 

Latent  heat  flux 

363  W  nr2 

436  W  nr2 

Buoyant  heat  flux 

275  W  m*2 

158  W  m'2 

Scaling  temperature 

-0.9664 

-0.5078 

Virtual  scaling  temperature 

-1.07 

-0.635 

Obukhov  length 

-3.52  m 

-5.63  m 

Scaling  ratio 

-0.5683 

-0.3551 

The  simulation  of  the  table  was  performed  by  using  a  surface  energy  balance  model  developed 
by  Rachele  and  Tunick  (1991). 

The  moisture  availability  coefficient  illustrated  in  figure  2,  determined  from  experimental  data, 
appears  to  be  a  valid  approach  for  weighing  the  potential  evaporation  as  a  means  for  estimating 
actual  evaporation.  These  preliminary  results  need  to  be  expanded  to  include  correlations 
between  daily  normalized  evaporation  rates  and  the  mean  daily  volumetric  water  content  of  the 
soil.  Other  verification  techniques  involving  the  obvious  relationships  between  evaporation 
versus  season,  soil  type,  latitude  and  climatic  extremes  are  obviously  needed. 

The  determination  of  the  latent  heat  flux  using  equation  (12)  yields  values  equivalent  to  the 
potential  latent  heat  flux  (equation  (14)),  indicating  that  the  moisture  availability  coefficient 
(equation  (16))  must  be  added  to  equation  (12),  such  that 


L  E  =  -  M  L  p  ut  q„ 


(23) 


resulting  in  an  expression  for  the  actual  latent  heat  flux. 
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5.  Conclusions 

The  impact  of  evaporation  and  the  presence  of  water  vapor  in  the  atmosphere  have  an  almost 
arcane  effect  upon  the  observed  dynamics  of  the  surface  boundary  layer.  The  addition  of 
moisture  to  surface  boundary  layer  models  such  as  the  dynamic  similarity  theory  is  in  reality 
quite  simple,  especially  if  the  moisture  availability  coefficient  approach  is  used.  The  apparently 
almost  independent  nature  of  the  coefficient  in  its  normalized  form  is  considered  to  be  a  valuable 
asset  to  the  modeling  of  atmospheric  processes.  As  additional  experimental  data  becomes 
available,  further  exploration  and  exploitation  of  the  moisture  availability  coefficient  will  be 
undertaken. 

E  =  MEp  +  (1  -  M)  Ea.  (24) 


Ritchie  (1972)  found  the  soil  limiting  rate  to  be 

E ,  =  C  t-'*,  (25) 


where  t  is  time  in  days  and  C  varies  with  soil  type  and  season  of  the  year.  Seasonal 
dependencies  are  a  function  of  temperature  and  will  vary  by  about  a  factor  of  2  from  winter  to 
summer  as  shown  in  figure  1.  Moisture  availability  coefficients  extracted  from  Jackson  et  al. 
(1976)  are  given  in  figure  2. 

The  complexity  of  the  evaporation  process  requires  that  the  energy  limiting  and  soil  limiting 
phases  be  partitioned  according  to  the  starting  times  of  the  soil  limiting  fraction.  Thus  equation 
(17)  must  be  rewritten  as 

E  -  MEp  +  C  £  (Af,_,  -  Af j)  (/»-*  +  1 Y112,  (26) 


where  n  is  the  number  of  days  after  evaporation  began. 
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Figure  3.  A  comparison  of  equations  (21)  and  (22)  in  terms  of  the  scaling  ratio  Z/L. 
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